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A pulse shape analysis algorithm and a method for suppressing the noise component of signals 
from a large copper proportional counter in the experiment aimed at searching for 2K capture of 
78 Kr are described. These signals correspond to a compound event with different numbers of charge 
clusters due to from primary ionization is formed by these signals. A technique for separating 
single- and multipoint events and determining the charge in individual clusters is presented. Using 
the Daubechies wavelets in multiresolutional signal analysis, it is possible to increase the sensitivity 
and the resolution in extraction of multipoint events in the detector by a factor of 3 -r- 4. 
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I. INTRODUCTION 

The main modern achievements in studying processes 
of double beta decay (/3~/3 - -decay) can be attributed 
to detection of its two-neutrino mode. This process has 
been discovered in as many as ten nuclei (see reviews 
[H j @ > El ) ■ The data obtained for the two- neutrino mode 
offer a chance to directly compare different models of the 
nuclear structure, which form the basis for calculations 
of nuclear matrix elements ||M 2l/ ||, and to select the opti- 
mal one. Though direct correlation between the values of 
nuclear matrix elements for the two-neutrino and neutri- 
noless modes of /3j3 decay is absent lacking, the methods 
for calculating ||M 2l/ |j and ||M 0ly || are very close, and a 
chance possibility to estimate their accuracy in calculat- 
ing ||M 0,y || appears only when comparing experimental 
data and theoretical results calculations for the proba- 
bility of 2fi(2i>) decay. 

It can be expected that acquisition of experimental 
data on the other types of 2/3 transitions (2/3 + ,K/3 + , 
and 2K processes) will make it possible to considerably 
increase the quality of calculations for both 2v and Ov 
decays. Much efforts have been currently made in search- 
ing for these processes ( [4] , 5] , [6] ) in spite of the fact that 
the 2[3 + (2v) and Kj3 + {2v) modes are strongly suppressed 
relative to 2/3(2^) decay due to the Coulomb barrier for 
positrons, and a substantially lower kinetic energy at- 
tainable in such transitions. Positrons are absent in the 
final state of the 2K(2v) decay, and the kinetic energy 
of the transition may be rather high (up to 2.8 MeV), 
which dictates determines an increased probability of a 
decay. However, this process is also difficult to detect, 
since it is only characteristic radiation that is detectable 
in it. The state-of-the-art experimental limit on the 78 Kr 
half-life period with respect to the 2K(2v) capture is 
Ti/2 < 1.5 x 10 21 yr (90% C.L.) [%. The theoretical cal- 
culations based on different models predict the following 
78 Kr half-lives for this process: 3.7 x 10 21 yr Q, 4.7 x 10 22 
yr @, and 7.9 x 10 23 yr The last two values were 

obtained from the estimates of the 78 Kr half-life with re- 



spect to the total number of 2e(2^) captures including 
in view of the 78.6% fraction of 2K(2v) capture events 



which makes 78.6% ll|. From comparison of the exper- 
imental and theoretical results, it is apparent that the 
sensitivity of measurements has reached the lower limit 
of theoretical estimates. 



II. THE TECHNIQUE OF THE EXPERIMENT 

The 78 Kr(2e K , 2^) 78 Se reaction produces a 78 Se** atom 
with two vacancies in its if-shell. The technique for seek- 
ing this reaction is based on the assumption that the 
values of energies of characteristic photons, and of the 
probability that they will be emitted when the double 
vacancy is being filled, coincide with the corresponding 
values when two separate single vacancies in the K shells 
of two isolated singly ionized Se* atoms are being filled. 
In this case, the total measured energy is 2K a \, = 25.3 
keV, where K a b is the binding energy of a if electron in a 
Se atom (12.65 keV). The fluorescence yield upon filling 
of a single vacancy in the if-shell of Se is 0.596. The en- 
ergies and relative intensities of the characteristic lines in 
the K series are K = 11.22 keV (100%), K = 11.18 keV 
(52%), K = 12.49 keV (21%), and K = 12.65 keV (1%) 
hj. There are three possible ways for deexcitation of a 
doubly ionized K-sheW: 1) emission of Auger electrons 
only (e a , e a ), 2) emission of a single characteristic quan- 
tum and an Auger electron {K, e a ), and 3) emission of 
two characteristic quanta and low-energy Auger electrons 
(K,K,e a ) , with probabilities pi = 0.163, p2 = 0.482, 
and p3 = 0.355, respectively. A characteristic quantum 
can travel a long enough distance in a gas medium be- 
tween the points of its production and absorption. For 
example, 10% of characteristic quanta with energies of 
11.2 a nd 12.5 keV is absorbed in krypton at a pres- 
sure of 4.35 atm (p — 0.0164 g/cm 3 ) on a path 1.83 and 
2.42 mm long, respectively (the values of absorption fac- 
tors are taken from [Hj|). The paths of photoelectrons 
with the same energies are 0.37 and 0.44 mm, respec- 



2 



tively. They produce almost pointwise charge clusters of 
primary ionization in the gas. In case of the event with 
the escape of two characteristic quanta absorbed in the 
working gas and a single Auger electron, the energy will 
be distributed among three pointwise charge clusters. It 
is these three-point (or three-cluster) events possessing a 
unique set of features that were the subject of the search 
in[|. 

A large proportional counter (IPC) with a casing made 
of Ml-grade copper is used to detect the above consid- 
ered processes. The LPC has a cylindrical shape with 
inner and outer diameters of 140 and 150 mm, respec- 
tively; its section along the axis is schematically shown 
in FigJTJ A gold-plated tungsten wire of 10 fim in diame- 
ter goes along the LPC axis and serves as the anode. The 
potential of +2400 is applied to the wire, and the casing 
(the cathode) is grounded. Both ends of the anode are 
lead to the appropriate end cap flanges via high-voltage 
pressure-sealed bushings-ceramic insulators with a cen- 
tral electrode taken from spark plugs. 




71Mmm 



FIG. 1: Schematic view of the LPC in section along the 
anode wire: 1 - wire (collecting electrode), 2 - load-carrying 
insulator, 3 - cathode, 4 - tubular bulge of the anode. 

To reduce the influence of edge effects on the oper- 
ating characteristics of the counter, the end segments of 
the wire are passed through the copper tubes with dimen- 
sions of 3.0 x 38.5 mm, which are electrically connected 
to the anode. Gas amplification is absent on these seg- 
ments, and charges are collected in an ionization mode. 
With the fluoroplastic insulator, the distance from the 
working region to the flange is 70 mm. 

The working part of the LPC is 595 mm in length (the 
distance between the end caps of the tubes); therefore, 
the LPC s operating volume is 9,159 1. The total capaci- 
tance of the counter and the outlet insulator is 30.6 pF. 
The total resistance of the anode and two output elec- 
trodes is 600 Om. All detachable joints are sealed with 
indium wire. All nipple joints are sealed with fluoroplas- 
tic gaskets. The internal insulators are made of fluoro- 
plastic. Their thickness was selected so as to be the small- 
est possible in order to improve the degassing conditions 
during vacuum treatment of the counter and stabilize its 



operating characteristics in the course of measurements. 

The LPC is filled with a pure Kr sample to a total pres- 
sure of 4.51 att; no quenching or accelerating gases are 
added. Prior to filling, Kr is purified of electronegative 
impurities in a Ni/SiO*2 reactor. 

The LPS's signals are read out by a charge-sensitive 
amplifier (CSA) from one end of the anode wire. The 
CSA parameters have been selected so that the signal is 
transmitted with minimum distortions, and information 
of the spatial distribution of primary-ionization charges 
in a projection onto a counter radius is fully represented 
by the pulse shape. When amplified in an auxiliary am- 
plifier, the pulses arrive at the input of the digital os- 
cilloscope LA-n20-12PCI, the output data of which (the 
pulse waveform digitized with a frequency of 6.25 MHz) 
are recorded with a personal computer. The length of 
the scanning frame with a resolution of 160 ns is 1024 
points (163.8 fxs), of which 50 fxs is the "prehistory" and 
114 fj,s is the "history". 

The counter is calibrated through the wall of its cas- 
ing by 7 rays of a 109 Cd source [E~ — 88 keV; rela- 
tive yield 0.036 photons/decay). Figure [5] presents (0) 
the total pulse amplitude spectrum and the energy spec- 
tra of (1) single-, (2) two-, and (3) three-point events 
from the source located in the middle of the LPC length. 
The procedure for obtaining them from digitized pulses 
is described in what follows. The following factors make 
their contribution to the low-energy part of the spec- 
trum: characteristic radiation AgK afl (E s» 22 keV) from 
this source, which "survived" after passing through a 5- 
mm-thick copper wall; scattered radiation from the wall, 
which is in equilibrium with the characteristic radiation, 
and Compton electrons from scattering of 88-keV pho- 
tons in the gas with the escape of a Compton photon 
beyond the counter. 
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FIG. 2: Pulse amplitude spectrum from the external 109 Cd 
source located in the middle of the LPC length: (0) all events, 
(1) single-point events, (2) two-point events, (3) three-point 
events, and (4) escape peak of characteristic 7-ray photons. 

The 88-keV peak is wider on the low-energy side due 



to the contribution of 88-keV 7 rays scattered from the 
wall. The energy resolution of this peak, determined by 
its right half, is 6.5%. Peak 4 at an energy of 75.4 keV 
corresponds to the escape of Kr characteristic radiation 
{EKr a = 12.6 kcV) beyond the counter. 

The 88-keV full-energy peak contains events with dif- 
ferent internal structures. Quanta with this energy are 
absorbed in Kr mostly by photoeffect in the if -shell 
(86.7%). The photoeffect in other shells makes 13.3% 
(l3| . }14| . Filling of the vacancy in the K-sheW of Kr 
is accompanied by emission of characteristic radiation in 
66.0% of cases and Auger electrons in 34.0% The 
theoretical efficiency of characteristic radiation absorp- 
tion in the counter s working volume is 86.9%. There- 
fore, the photoeffect is responsible for 49.7% of two-point 
events (0.867 x 0.660 x 0.869) and 42.8% of single-point 
events (0.133 + 0.867 x 0.340) out of the total num- 
ber of absorptions due to photoeffect in the full-energy 
peak. By single-point events, we mean all events in 
which only electrons escaping from the shell of a sin- 
gle atom (photoelectrons -I- Auger electrons), including 
events of photoelectric absorption in the upper shells 
of Kr. Only single-point events in an amount of 7.5% 
[0.867 x 0.660 x (1.000 - 0.869)] of the total number of 
photoabsorption events will be presented in peak 4- 

Some primary quanta can be absorbed as a result of 
two-step process of " Compton scattering-photocffect" . A 
Compton electron creates one ionization point. A Comp- 
ton photon absorbed by photoeffect participates in the 
above-described processes. Therefore, the two-step pro- 
cess makes its contributions to the full-energy peak in 
the form of two- and three-point events and to peak 4 in 
the form of two-point events. Upon normalization to the 
peak area, the estimated final composition of events for 
the full-energy peak contains 44.1% (single-point events) 
+ 51.2% (two-point events) due to photoeffect + 2.2% 
(two-point events) + 2.5% (three-point events) due to 
the two-step process. In peak ^, there are 95.3% single- 
point events + 4.7% two-point events. 

Examples of pulses (dark lines) corresponding to events 
of two types are presented in Fig[3l The pulse due to pho- 
toabsorption of a 88-keV photon with escape of electrons 
only (a single-point event) is shown in Fig|3^, and simul- 
taneous escape of characteristic photon Kr with an en- 
ergy of 12.6 keV and a photoelectron (E 1 — E) = 88 keV 
- 12.6 keV=75.4 keV (a two-point event) is illustrated 
in Figj3j). The maximum distance between pointwise 
charge clusters in projection onto the counter s radius 
is equal to the radius. For pure Kr, the evaluated time 
it takes for ionization electrons to drift from the cathode 
to the anode is 53 fj,s. From Figsj3^i and[3jD, it is appar- 
ent that the second pulse with a ~ 3 4- 5 times smaller 
amplitude is produced in the counter in about ~ 53 /j,s 
after the first one. This pulse is generated by secondary 
photoelectrons knocked out of the cathode by photons 
produced during development of an avalanche from the 
primary ionization. The probability of photoeffect on the 
cathode is rather high, since the working gas contains no 
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FIG. 3: Examples of pulses (dark lines) of two types of events: 
(a) for photoabsorption of an 88-keV photon with an escape of 
electrons only (a single-point event), (6) simultaneous escape 
of a 12.6-keV characteristic photon Kr Ka and a photoelectron 
(E 7 - E A >A' a )=88.0 keV-12.6 keV = 75.4 keV (a two-point 
event). The calculated, area-normalized current pulses of pri- 
mary ionization electrons are shown with light lines in Figs. 
3c and 3d, while the respective voltage (charge) pulses ob- 
tained by integrating these current pulses are depicted with 
light lines in Figs. 3 a and 36. 



quenchers. 



III. DETERMINING THE SHAPE OF THE LPC 
CURRENT SIGNAL 

The digitized output pulse from the measuring channel 
can be presented as direct convolution of sought signal 
x(U) with response pulse G from the linear system, which 
is distorted by stochastic or determinate noise z(U): 



y(ti) := Gx(ti)+z(ti),i = 0,...,N-l 
:= (g 03) x)(ti) + z{U). 



(1) 



The noise level determines the lower limit on the sen- 
sitivity of the measuring channel, while the instrumental 
function defines the resolution value. 

Given the values of y and G, one can try to estimate 
x(t) in the presence of noise z(t); i.e., in principle, it is 
possible to state the inverse problem (deconvolution) of 
determining the signal at the linear system output by the 
values of the output signal: 



x(ti) = G^yiU) = x{U) + G~ 1 z(t i ), 



(2) 



where G _1 is the operator inverse. 

To determine the values of charges released in indi- 
vidual clusters of a multipoint event, one can differenti- 
ate the original charge pulse ^G _1 — -§^j and represent 

the obtained shape by a set of Gaussian curves. The 
evaluated area under an individual Gaussian curve will 
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correspond to the charge (energy) value in the relevant 
cluster. From Figs[3]c and[3]i, it is apparent that direct 
differentiation provides an asymmetric bell shape (dark 
line) . Such shape results from the nearly Gaussian distri- 
bution of the current pulses due to electrons of primary 
ionization from a pointwise energy deposit, which arrive 
at the boundary of the gas amplification region near the 
anode wire. 

This shape is determined by the spatial distribution of 
the charge density in projection onto the radius. The pa- 
rameters of this distribution depend on the time it takes 
for the primary charge cluster to drift to the anode. As 
it drifts, the charge cluster spreads out into a cloud due 
to electron diffusion. The pulse read out from the anode 
wire is mostly produced by a negative charge induced on 
the anode moving toward the cathode by positive ions 
produced near the wire in gas amplification process and 
moving toward the cathode i.e., the ion component (i.e.). 
The estimated total ion drift time is 0.447 s. The contri- 
bution of the equilibrium (with ions) electron component 
(e.c.) to the total induced charge is ~ 7%. The electron 
collection time is ~ 1 ns. 

The output pulse shape is defined by the superposition 
of induced charges from single electron avalanches dis- 
tributed in time and intensity according to: the shape of 
the current pulse from primary ionization electrons, the 
shape of a pulse from an individual avalanche, and a finite 
time of the CSA self-discharge. The last two parameters 
are responsible for the asymmetry of the output current 
pulse. The output current pulse can be transformed to 
a symmetric shape by taking into account the analytical 
dependence of the amplitude of the output voltage pulse 
generated by a point (in projection onto the radius at 
the boundary of the gas amplification region) group of 
primary ionization electrons as a function of time and 
discharge constant of the output storage capacitor : 

V k (ti) = Ki. c n k cxp f~ ~^~r~ ^ x ( 3 ) 

X {^( 1 + l) + ife + 2^|^ + --j + 
+K e . c .n k exp (-j^jj i 

where Vk(U) is the amplitude of the voltage pulse from 
the fcth group of electrons, n& is the number of primary 
electrons in the kth group, U = t(to + tio) is the current 
time for the dependence of the voltage pulse amplitude 
from the kth group, to is the time of origin of the total 
pulse, and Uq is the time of origin of the pulse from the 
fcth group, ifi.c. = MV(i)j. c ,, M is the gas amplification 
factor, 

V(i)i.c. = (e/C) ■ln(r k /r )\n(r k /r a ) 

is the total voltage pulse amplitude produced at output 
capacitor C by a single ion generated in the gas dis- 
charge, e is the electron charge, r$ is the radius corre- 
sponding to the avalanche's center of gravity, r a is the 



anode radius, r\ is the cathode radius, B is the time 
parameter associated with the motion of positive ions 
of the gas discharge in a particular gas (for the LPC 
filled with Kr at 4.51 at, B = 2.28 ns), RC = r d is 
the CSA discharge constant, R is the leakage resistance, 

K c . c . = MV"(i) e .c, and 

V(l)e.c. = (e/C) • ln(ro/r a )ln(rk/r a ) 

is the total voltage pulse amplitude produced at output 
capacitor C by a single electron generated in the gas 
discharge. In Eq.(3), the electron component is assumed 
to appear instantly. 

At Td — oo, Eq. (3) assumes the form 

V k {U) = Ki c ,n k x 'nf 1 + ^j+ K c . c .n k . (4) 

In our case, Td ra 192 fis. For a time interval satisfying 
condition t/rj < 2, Eq. (3) can be reduced to the first 
two expansion terms. 

If its is assumed that the gas discharge from nu of pri- 
mary electrons happens at the beginning of the digitiza- 
tion interval, the pulse amplitude at the end of this inter- 
val can be described by Eq. (3) for t — 160 ns, since the 
influence of the output capacitor discharge over this time 
is negligible. If a contribution of the earlier discharges is 
absent in this time interval, the pulse amplitude at the 
upper bound of the interval can be used to determine 
the rife value. In this case, it is taken into account that, 
at the end of a 160-ns interval, the contributions of the 
terms to the total Vk(U) value in Eq. (3) make 74 and 
26%, respectively. These conditions are satisfied in the 
recorded actual pulse in the first time channel from the 
beginning of the pulse. The n\ value obtained from the 
actual pulse is used in Eq. (3) to calculate the total shape 
of the partial pulse in the entire time interval from the 
beginning to the end of the frame. The pulse obtained 
thereby is subtracted from the actual one. Therefore, the 
above condition is now fulfilled for the first digitization 
interval of the residual pulse or for the second interval of 
the original pulse. This procedure is repeated until the 
last time channel in the frame. 

The sequence of nk values for a single-point event has a 
symmetrical distribution with a nearly Gaussian shape. 
It is this distribution that is used for further analysis. 
The area under the Gaussian curve or, in the case of a 
multiparticle event, the sum of the Gaussian areas on a 
time interval of 53 /is from the beginning of the pulse 
yields the total number of primary ionization electrons. 
To plot the spectra in Fig. 2, this sum is multiplied by 
the coefficient equal to the averaged ratio of areas of the 
actual current pulse and of the calculated Gaussian curve 
for purely single-point events. 

The calculated, area-normalized current pulses of pri- 
mary ionization electrons at the boundary of the gas am- 
plification region are shown with light curves in Figs[3]c 
and [3]ff, and the corresponding voltage (charge) pulses 
obtained by integrating these current pulses are depicted 
with light lines in Figs and 



5 



From Figs{3]c and [3]cZ, it is apparent that, at an en- 
ergy deposit of 88 keV, the signal-to-noise ratio (SNR) is 
rather high. In the energy range of 20 -j- 30 keV, in which 
the 2if-capture in 78 Kr is sought, the SNR for individual 
components of the total energy deposit corresponding to 
the possible effect (25.3 keV) is not so favorable. 

FiguresHJj andS^ present examples of recalculated cur- 
rent pulses for two types of two-point energy deposits due 
to the if-capture of 81 Kr isotope. The total energy de- 
posit corresponds to the binding energy of an electron 
in the if -shell of a daughter 81 Br (13.5 keV). The ener- 
gies of a characteristic quantum (Eku = 11.9 keV) and 
concomitant Auger electrons (E a = 1.6 keV) are close to 
the energies of individual components for events of the 
2i4f-capture of 78 Kr. More detailed information on the 
81 Kr source is given in what follows. 




FIG. 4: Primary ionization signals for two types of two-point 
energy deposits from Tf-capture of 81 Kr isotope accompanied 
by an escape of a characteristic photon and a concomitant 
Auger electron: (a, b) recalculated voltage (charge) pulses and 
(c, d) corresponding current pulses from primary ionization 
electrons. 

From Fig 01 it is apparent that these signals have a 
high noise level. Noises and possible electric pickup may 
both mask the low-energy component and create a false 
one. The use of traditional methods of frequency filtering 
with different window functions, e.g., in the form of the 
Hamming [l6j], Wiener pj}, and Savitzky-Golay [HI fil- 
ters, sometimes fail to ensure reliable extraction of closely 
spaced (Fig. 4a) and masking each other components of 
a compound event. 



IV. 



WAVELET-BASED FILTER 



Mathematical studies carried out in late 1980s initiated 
intense development of a principally new class of orthog- 
onal transforms based on the use of wavelet functions 
19 . Wavelet transforms are distinguished by a high de- 
gree of locality of their base functions both in the time 
and frequency regions, which allows one to use them for 



processing of many nonstationary processes. "In the pre- 
liminary processing of our LPC data it is expedient to use 
multiresolutional signal analysis [2(| based on the dyadic 
transform of discrete signals, and often called discrete 
wavelet transform (DWT). In this case, the analyzed sig- 
nal x(t) is presented as the decomposition 



X (t) — ^ ] J j .m l Pjo,m 
m=0 



(*)■ 



j=j m=0 



i>j,m(t), (5) 



where are the well-known orthonormal scaling (scaling 
functions) and wavelet functions- "ancestor" wavelets; 
Ojo,m — ^iVjo,m) ar e the empirical approximation co- 
efficients and dj m — (x, ^j,m) are the empirical detailing 
coefficients; j,m € Z are the current values of the scale 
and the shift; Nj (Nj) is the number of approximation 
(detailing) coefficients considered at the relevant levels of 
the decomposition; jo is the initial scale value; and J is 
the final scale value. 

Parameter J specifies the resolution of wavelet recon- 
struction x J (t) of original signal x(t). Actually, in the 
case of tendency J —> oo, the norm tends to \\x J — x\\ —> 0. 

Scaling functions {4>j. m (t)} and mother wavelet func- 
tions {4>j, m {t)}j,mez, which have 2K nonzero coefficients, 
satisfy the so-called two-level relations [ll| : 

2K-1 

tp(t) = a/2 ^2 h(n)tp{2t - n), 

n=Q 



2K-1 

i/>(t) = V2 9{n)y{2t - n), 

n=0 

where h(n) and g(n) are the coefficients of low- and 
high-frequency filters of the wavelet transform, when 
g(n) = (— l) n h(2K — n — 1). As distinct from other 
types of transforms, in which the base functions are ex- 
plicitly specified, one succeeds in analytically obtaining 
the base functions in the wavelet analysis only in rare 
cases, and the basis is most frequently specified by co- 
efficients h(n) and g(n). In this paper, we use wavelets 
of the Daubechies family. The scaling functions and the 
Daubechies wavelets are continuous functions that are 
not identically equal to zero on a finite segment and are 
not differentiable anywhere on this segment. Table U con- 
tains the filter coefficients used in the Daubechies scaling 
functions db4 and db6 (numbers 4 and 6 denote the num- 
ber of nonzero coefficients in the filters). They are ra- 
tional numbers and fully define the Daubechies wavelet 
transform (DVT) [l9[. In contrast both to the Fourier 
transform, which localizes frequencies, but fails to pro- 
vide time resolution for the process, and to the apparatus 
of (5-functions, which localizes moments of time, but have 
no frequency resolution, the DVT makes it possible to re- 
veal the local properties of any structure of an individual 
event on different scales, eliminating smooth polynomial 
characteristics and emphasizing fluctuation structures. It 
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TABLE I: Coefficients of the low-frequency filters for wavelets 
with compact carriers db4 and dbQ. The h(n) values are nor- 
malized so that Xln=o 1 M n ) = V^2 (191 ] 



h(n) 



h(0) 
h(l) 
h(2) 
ft(3) 
h(4) 
h(5) 



db4 (K = 2) 



(3 + V3) /4\/2 
(3 - V3) / 4>/2 
(l-\/3)/V2 



rib6 (X = 3) 



1 + VW + V5 + 2VW) / lQy/2 



5 + \/T0 + 3^5 + 2^/10) / 16^2 



10 - 2^10 + 2\/5 + 2yT0) / 16^2 



10 - 2^/TO - 2V5 + 2^/10 J / 16^2 
5 + VlO - 3\/5 + 2vTo) / lGv^ 
1 + \/l0 - y/5 + 27Io] / 16^2 



is possible to completely eliminate statistical fluctuations 
by selecting only strong correlated fluctuations, which 
will allow one to observe exactly those dynamic fluctua- 
tions that exceed the statistical component (noise). 



A. Noise elimination using the wavelet method 

Among the modern noise suppression techniques, the 
Donoho- Johnstone method has received the widest ac- 
ceptance . This method is rather simple in implemen- 
tation and time-saving in computational aspect, since it 
implies the use of only fast algorithms of wavelet trans- 
form. It consists of three steps, which, being successively 
applied to the original signal, produce a noise-suppression 
effect. At the first step, signal under investigation x(t) is 
subjected to DWT; afterward, the threshold noise elimi- 
nation procedure is applied to each of the detailing coef- 
ficients of level j and, sometimes, to the approximation 
coefficients of the same level; and, finally, inverse wavelet 
transform is performed, which results in reconstruction of 
the signal that is characterized, as expected, by a higher 
SNR. 

This technique is a nonparametric estimate of the re- 
gression signal model with the use of an orthogonal ba- 
sis [12] ~ [3]; it operates with high efficiency with signals 
in the decomposition of which only a few detailing co- 
efficients significantly differs from zero. Selection of a 
particular wavelet form depends on the problem under 
investigation and is not predetermined beforehand pBI |. 
The depth of decomposition generally depends on the 
properties of the analyzed signal. Smooth wavelets pro- 
duce a smoother signal approximation, and vice versa- 
" short" wavelets better search for peaks of the approxi- 
mated function. The depth of decomposition affects the 
scale of rejected details; i.e., as the depth of decomposi- 
tion increases, the model subtracts the noise of a steadily 
increasing level until the scale of details becomes too 
large and the transform starts distorting the original sig- 
nal shape. Upon further increase in the depth of de- 
composition, the transform starts forming a smoothed 
version of the original signal; i.e., apart from the noise, 
some local peculiarities of the original signal are also fil- 



tered off. Decomposition of the signal and its reconstruc- 
tion by the approximation and detaining coefficients that 
have passed threshold processing is carried out using the 
Malla algorithm [H and the lifting procedure [Hj]. The 
threshold processing itself is performed with the aid of 
one of the threshold noise suppression operations. In ac- 
cordance with the hard threshold processing [2lj . all co- 
efficients \k € z \ of level j, which are greater than 

or equal to the threshold, are held constant, while the 
other coefficients that do not satisfy this condition are 
are reduced to zero: 



4* = ( 



ij,k 



> 



(6) 



where pj is the threshold factor for the specified scale [20[ , 
and <t is the noise variance on the jth scale. Hard thresh- 
old processing is a straight-out procedure. It should nev- 
ertheless be taken into account that a hard threshold 
has two drawbacks that can lower its usefulness for the 
noise suppression task. The first consists in the fact that 
retaining of detailing coefficients above a predetermined 
threshold value also implies retaining of their noise. The 
other drawback is in the presence of parasitic harmonics 
generated in the resulting signal by artificial introduction 
of lacunas (gaps) formed due to the coefficients reduced 
to zero. The use of soft threshold processing implies recal- 
culation of detaining coefficients d in the following way: 



sign (dj,kj ( 



Jtrh 



(7) 



where is a certain threshold value. 

In this case, apart from the reduction to zero of the co- 
efficients d actually containing only the noise component, 
the detailing coefficients are decreased by the Othr value, 
which corresponds to noise suppression in the informa- 
tive coefficients as well. Threshold Othr, which defines 
the analyzed signal, is a sole parameter that cannot be 
exactly estimated directly from actual data. Mathemat- 
ical solutions in determining the threshold value seem 
unsatisfactory to be a priori accepted. In an ideal case, 
we tend to specify this threshold so as to minimize the 
root-mean-square error (rmse) : 



\ 



1 N 

- b(*0 - x{u)Y 



(8) 



Therefore, by analogy to [20|,[26|, we determine the ba- 
sic signal component relying on threshold processing of 
the model signal. Ideal threshold processing minimizes 
the rmse, making all coefficients with the signal compo- 
nent below a vanish [20l |. In other words, it is assumed 
that the free from noise signal component above a is an 
original signal; i.e., Othr = cr. This selection can be tested 
on the model signals. Our experience gained in working 
with natural signals confirms that the ratio of the sig- 
nal peaks to the noise level reaches its maximum near 
Othr = cr. 



B. Test of model signal with noise 

To illustrate the efficiency of different approaches to 
the deconvolution problem, let us select model signal 
x(t n ), which is specified on a set of points N — 2 10 and 
consists of three closely spaced Gaussian peaks 

a k exp (-(t - t k ) 2 /2fil) , 

where a k = [1.32,7.0,7.0], t k = [277,307,327], and fi a = 
7. Let us add Gaussian noise and a harmonic term to it: 

x(t n ) = x(t n ) + kA(t) + g(t n ), 

where A(r) is ordinary Brownian process analyzed on r £ 
[0, 1], k is a certain coefficient, and is the harmonic 
term. 

For the variance <r 2 of the noise component in model 
signal W(t n ) — kA(t) + g(t n ) to be estimated by the 
x{t n ) data, it is necessary that the influence of signal 
x{t n ) be suppressed. A rough estimate can be obtained 
from the mean values of the smallest-scale wavelet coef- 
ficients [2l[. The signal with length N has N/2 wavelet 
coefficients {(x(t n ), "4>j.m)}Q <m<N / 2 °^ * ne smallest scale 
2 J = 27V~ 1 . Coefficients |(x(t„),-0j,m)| are small if sig- 
nal x(t n ) is smooth on carrier ipj t m, and in this case, 
(x(t n ), ipj >m ) PS (W,ipj >m ). However, if x{t n ) has a sharp 
differential on carrier ip, coefficients \(x(t n ),ipj,m)\ are 
rather large. In the case of a piecewise-smooth signal, we 
have several differentials, which results in a certain num- 
ber of large coefficients. This number is small in compar- 
ison with N/2. At the smallest scale, signal x(t n ) defines 
the value of a small portion of large-amplitude coefficients 
\{x(t n ),ipj tm )\, which are considered to be "overshoots". 
All the other coefficients are approximately equal to 
(W,ipj im ); they are independent random Gaussian vari- 
ables with variance a 1 . Therefore, a 1 can be roughly 
estimated, with the influence of x(t n ) being ignored, by 
the mean of absolute values {(x(t n ), V'j,m}}o< m <jv/2' 
vided by 0.6745 [H[. 

Daubechies filters db4 were used in a direct wavelet 
transform and when the estimate of the useful signal was 
synthesized from coefficients after their threshold pro- 
cessing, whereas higher-order filters db6 were employed 
to estimate the threshold values. 

Calculations of all wavelet coefficients, testing of the 
noise suppression algorithms, and selection of required 
signals for further analysis were performed using MAT- 
LAB and its web libraries 271 ] . 

The results are shown in FigJS] It is apparent that 
all filters provide good noise suppression efficiency (the 
SNR is > 30 dB), but the signals passed through the 
wavelet filters with the use of the hard or soft threshold 
processing better represents the shape of the original sig- 
nal (FigsJSJs and[5|f), which ensures more reliable signal 
resolution into individual components. For our purposes, 
hard threshold processing is quite satisfactory both in the 
quality factor and in the computational time. 
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FIG. 5: (a) Model signal from three closely spaced Gaussian 
peaks; (b) additive mixture of the model signal with noise; 
noise suppression using (c) the Savitzky-Golay filter and (d) 
Wiener filter; (e, /) Daubechies wavelet filters, db4, with hard 
and soft threshold processing, respectively. 



RELIABILITY IN IDENTIFYING 
MULTIPOINT EVENTS 



In order to verify the efficiency of the algorithm for 
the noise suppression and pulse shape discrimination for 
multipoint events from an actual source the LPC was 
filled with natural Kr containing cosmogeneous radioac- 
tive isotope 81 Kr (Ti/ 2 = 2.1 x 10 5 yr) with a bulk ac- 
tivity of - 0.1 min-H- l Kx [Hj], [H. It decays by elec- 
tron capture into 81 Br*. In 87.5% of cases, an electron 
is captured from the iGshell of a Kr atom (if-capturc) 
30] . Filling of the vacancy at the iGshell of a daughter 
Br atom is accompanied in 61.4% of cases by an escape 
of two characteristic quanta with energies of 11.92 keV 
(K al , 100%), 11.88 keV (K a2 , 50.9%)M3.29 keV (JG31, 
21.0%), and 13.47 keV (Kp 2 , 1.07%) \V2} and concomi- 
tant Auger electrons with energies of 1.55 keV, 1.60 keV, 
0.27 keV, and 0.01 keV, respectively (the relative inten- 
sities of K aj p lines are presented in brackets). From this 
list, it is apparent that, when photons K a \ and K a 2 es- 
cape, the energy of Auger electrons is sufficient for a dis- 
tinguishable two-point event to be formed. In 38.6% of 
cases, filling of the vacancy at the iGshell of Br is accom- 
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FIG. 6: Amplitude spectra (1) of single-, (2) two-point, and 
(3) three-point events that form the total analyzed pulse am- 
plitude spectrum of the counter detecting both background 
and 81 Kr decays. 



panied by an escape of a cascade of Auger electrons that 
produce a single-point energy deposit. Events with an es- 
cape of characteristic photons Kpi and K@2 should also 
be placed into a category of single-point events. Taking 
into account the absorption efficiency for characteristic 
photons in the LPC working gas (e a bs — 0.869) and us- 
ing the above data, one can calculate the composition 
of a 13.5-keV full-energy peak: 49.4% of all events are 
single-point, and 50.6% are two-point. 

Describing the noise-free signals by a set of Gaussian 
curves using rmse minimization technique (7), one can 
discriminate between two- and single-point events. Fig- 
ure [6] presents the energy spectra (1) of single-, (2) two- 
, and (5) three-point events forming the total analyzed 
pulse-amplitude spectrum of the counter that detects 
both background and 81 Kr decays. Similar components 
(spectra 1-3) of total spectrum 1 for a calibration mea- 
surement of 88-keV 7 rays are shown in Fig. 2. 

Given the composition of events from the 81 Kr and 
109 Cd sources and results of simulation of these processes 
in the LPC, one can, using the GEANT 4.8.2 program 
[3l| . determine the quality of suppression of the noise 
component and the efficiency of the procedure for sepa- 
rating events by the criteria of a multipoint event. Table 
II presents the percentage composition of three types of 
events in the full-energy peaks with energies of 13.5 keV 
( 81 Kr) and 88 keV ( 109 Cd), which was obtained by es- 
timation (column /); by simulation of processes in the 
LPC, ignoring the confluence of closely spaced primary 
charge clusters (column Ha), and taking into account 
their confluence (column lib); and by separation of the 
experimental spectrum into components after its wavelet 
purification of noise using hard threshold processing (col- 
umn III). 

Determining the parameters of individual components 
in a compound event, it is possible to perform energy 
calibration in the energy range under investigation and 
determine the energy resolution of single-point compo- 
nents which form a multipoint event. The energy dis- 



tributions of individual components of two-point events 
from the 81 Kr and 109 Cd sources are shown in Figs. 7 
and 8, respectively. 



100 

80 



% 60 



40 



(88-12.6) 



12.6 



1500 1000 500 



intensity 




1500 
I' 1000 
S 500 


20 40 60 80 100 
E , keV 

FIG. 7: Amplitude distributions of energy deposits of individ- 
ual components for two-point events from the source placed 
in the middle of the LPC length ( 109 Cd). 
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FIG. 8: Amplitude distributions of energy deposits in of indi- 
vidual components of for two-point events from the internal 
( 81 Kr) source with a total energy of 13.5 keV. 
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VI. CONCLUSIONS 

The analysis both of the pulse shape variety and meth- 
ods of LPC signal processing has been performed using 
a multiresolutional wavelet analysis. It has been shown 
that the noise component can be suppressed successfully, 
and the problem of separating of the events with closely 
spaced clasters, partially or fully lapping over each other, 
can be solved with a high efficiency, and the parameters 
of these events can be estimated. The test measurements 
of signals from the krypton-filled LPC irradiated by in- 



ternal and external sources have made it possible to de- 
termine the efficiency of discrimination between single-, 
two-, and three-point events. The described technique 
helps increase the sensitivity of detection of compound 
events in gas-filled detectors and demonstrates its appli- 
cability to detection of rare events from 2K{2v) capture 
in large proportional counters. 
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TABLE II: Relative composition of three types of events in the full-energy peaks with energies of 13.5 keV ( Kr) and 88 keV 
( 109 Cd), obtained by estimation (column 7); by simulation of processes in the LPC, ignoring the confluence of closely spaced 
pointwise ionization regions (column Ha) and taking into account their confluence (column lib); and by separation of the 
experimental spectrum into components after its wavelet purification of noise using hard threshold processing (column 777). 





13.5 keV ( sl Kr) 


88 keV ( I,JM Cd) 




calc, % 


GEANT4.8.2, % 


empiric, % 


calc, % 


GEANT4.8.2, % 


empiric, % 


type of events 


7 


77a 


lib 


III 


7 


77a 


lib 


III 


singl-point 


49.4 


44.1 


66.2 


68.2 


44.1 


45.0 


59.4 


56.8 


two-point 


50.6 


55.9 


33.8 


31.7 


53.4 


53.7 


40.0 


39.4 


three-point 








< 0.1 


2.5 


1.3 


0.6 


3.4 
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